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ABSTRACT 

Recently is has been suggested that the dearth of smaU halos around the Milky Way arises due to a 
modification of the primordial power spectrum of fluctuations from inflation. Such modiflcations would 
be expected to alter the formation of structure from bottom-up to top-down on scales near where the 
short-scale power has been suppressed. Using cosmological simulations we study the effects of such a 
modification of the initial power spectrum. While the halo multiplicity function depends primarily on 
the linear theory power spectrum, most other probes of power are more sensitive to the non-linear power 
spectrum. Collapse of large-scale structures as they go non-linear regenerates a "tail" in the power 
spectrum, masking small-scale modifications to the primordial power spectrum except at very high-z. 
Even the small-scale {k > 2 h Mpc~^) clustering of the Ly-a forest is affected by this process, so that 
CDM models with sufficient power suppression to reduce the number of lO^^M© halos by a factor of ~ 5 
give similar Ly-a forest power spectrum results. We conclude that other observations that depend more 
directly on the number density of collapsed objects, such as the number of damped Ly-a systems, or the 
redshift of reionization may provide the most sensitive tests of these models. 

Subject headings: cosmology: theory - large-scale structure of universe 



1. INTRODUCTION 

The field of physical cosmology has made rapid progress 
in the last decade, and a "standard model" is already be- 
ginning to emerge. Many of the main cosmological pa- 
rameters are becoming known and there is good reason 
to believe that the measurements will be significantly im- 
proved, and the paradigm tested, in the next few years 
from observations of the Cosmic Microwave Background 
anisotropy and upcoming surveys of large-scale structure. 
While in broad outline the paradigm appears to work well, 
there are some discrepancies which indicate that revisions 
in our standard model may be required. In this paper we 
discuss several topics related to one of these issues: the 
lack of low mass halos in our local neighborhood and in 
particular consider what we might learn about the small- 
scale matter power spectrum. 

The halo problem has been highlighted by sev eral gr oups. 

Analytic arguments based on Press-Schechter 

ory were given by Ka uffma nn, White & Guiderdoni (|1993|), 



while Klypin et al. ( |1999D and Moore et al. ( |l999aD used 
very high resolution dark matter simulations. A summary 
of the s ituat ion has been given recently by S perge l & Stein- 
hardt ( |1999D and Kamionkowski & Liddle (|1999| ). 

Within the Press-Schechter theory, and its extensions, 
the number density of halos of a given mass is related to 
the amplitude of the linear theory power spectrum on a 
scale proportional to M^^^. For example 10^°Mq halos in 
a model with = 0.3 probe a linear scale of 0.3ft.~^Mpc. 
Numerous numerical simulations have demonstrated that 
halo number density seems to be governed by the linear 
power spectrum as Press-Schechter theory would predict. 
A deficit of low mass halos thus implies cither additional 
physics (see below) or a deficit of linear theory power on 
small length scales. This modification could come about 
either by variations in the primordial power spectrum 
(e.g., from infiation) or in the cosmological processing of 
this power spectrum (e.g., from 'warm' dark matter). 



Recently Kamionkowski & Liddle ( 199£ ) pointed out 
that a well studi ed cla ss of infiationary models (BSI; see 
e.g., Starobinsky 1992) could give rise to a deficit of small- 
scale power in the primordial power spectrum. One way to 
achieve such a deficit is to introduce a change in the slope 
of the inflaton potential at a scale determined by the astro- 
physical problem to be solved, in our case k ~ 5/iMpc~^. 
Thus in such models one 'naturally' achieves fewer low 
mass halos than in the conventional inflationary CDM 
models. 

The linear theory power spectrum of these BSI mod- 
els is well described by Kamionkowski & Liddle ( 1999 ). 
Compared to a scale-invariant model there is a small rise 
followed by a sharp drop in power at some scale feg. Be- 
yond ko the power spectrum oscillates with an envelope 
which falls more steeply than k~^. At high-fc the spectrum 
recovers to the usual k~^ slope but with much smaller am- 
plitude. 

We believe that it is of interest to constrain such modifi- 
cations to the initial power spectrum, if possible. However, 
because of this sharp drop, the model near feg more closely 
resembles the familiar top-down scenarios (e.g. HDM) than 
a bottom-up CDM model for some range of wavenum- 
bers. Thus arguments based on reasoning developed for 
'traditional' CDM models should be checked against nu- 
merical simulations. Furthermore, the number density of 
objects may be one of the only probes of the linear theory 
power spectrum. Several astrophysical probes of small- 
scale power are sensitive not to the linear theory but to 
the non-linear power spectrum. As is well known, objects 
collapsing under gravitational instability feed power from 
large scales to small thus allowing small-scale power to 
be regenerated on ce a mode goes non-linear (e.g. Little 
Wein berg & Park |1991| ; Melo tt fc Shandarin |1990| , |1991 



1993; Bagla & Padmanabhan 1997). How much power is 



regenerated, crucial for determining how much there was 
initially, requires numerical calculation. 

We address several of these issues in the following sec- 
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tions. 

Finally we should note that we believe that constraints 
on small-scale power are of intrinsic interest in and of 
themselves. We shall discuss this within the context of 
the sub-halo problem described above while noting that 
several other astrophysical effects may also explain the 
discrepancy. The most obvious examples are: the to- 
tal number of Local Group satellites could be underesti- 
mated, feed back could be important (Kauffmann, White & 
Guiderdoni 1993, Bullock, Kravstov & Weinberg 2000| ), or 
the satellites could fail to make stars and be dark (e.g. HI 
clouds) . 

2. PROBES OF SMALL-SCALE POWER 

Any proposal to solve the small halo number density 
problem by modifying the initial power spectrum must 
simultaneously be able to pass other constraints on small- 
scale power. While we have a number of constraints on 
the linear and evolved power spectrum on larger scales, 
there are very few stringent constraint s on l inear scales 
below a Mpc. Kamionkowski & Liddle ( 1999| ) argue that 
constraints from the abundance of damped Ly-a systems 
and the reionization epoch are passed by low-density ver- 
sions of the BSI model, partly because of the uncertainties 
involved in making those predictions. The clustering of 
objects at high-z does not appear to be a promising probe 
of the matter power spectrum on these small scales. A pri- 
ori the two most obvious probes are the object abundances 
which motivated this modification of the power spectrum 
initially and the power spectrum of the flux in the Ly-a 
forest. 

3. SIMULATIONS 

To address some of these issues we ran two sets of 
N-body simulations. The base model in all cases was 
a ACDM model with a„ = 0.3, CIa = 0.7, h = 0.7, 
n = I, C OBE normalized using the method of Bunn 
& White (|1997|), i.e. with as = 0.88. The transfer 



func tions were computed using the fits of Eisenstein & 
Hu ( 199S| ) without the baryonic oscillations. This power 
spectrum was optionally filtered to suppress small-scale 
power. We have modeled the behavior displayed in Fig. 1 
of Kamionkowski & Liddle (1999) with a simple analytic 
form: 



A^(fc)^ 



_ k^P{k) 
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A^^ + ik/kof^'A^^iko))' (1) 



Here A"^^ is the fiducial power spectrum whose high-/c be- 
havior we are modifying and the power-law slope of the 
fc/fco term was chosen to match the behavior of their Fig. 1 
in the range just above fcg. The model plotted in their 
Fig. 1 corresponds to kg ~ 10/iMpc~^, as shown in our 
Fig. H below. 

The first set of simulations used a P M cod e descri bed i n 
detail in (Meiksin, White & Peacock |T999| , White |1999D . 
The simulations used 256^ particles and a 512'^ force mesh 
in a box 25/i~^Mpc on a side evolved from z = 70 to 
z = 3. The high mass resolution and quick execution times 
allowed us to explore parameter space and address the Ly- 
a forest questions (§ |4.3| ) where very high force resolution 
isn't necessary. 

The second set of simulations used a new implemen- 
tation of a TreePM code similar to that described in 



Bagla ( 1999 ). These runs used 128'^ particles in the 
same size box, evolved from z = 60 to z = 3 with 
the time step dynamically chosen as a small fraction of 
the local dynamical time. While higher mass resolution 
would be preferable, this would make the execution time 
prohibitive on desktop workstations with the current se- 
rial version of the code. A spline softened f orce (Mon- 
aghan & Lattanzio 1985, Hernquist & Katz 1989) with 



ft, = 8 X lO^^'Lbox — 20ft^^kpc comoving was used (the 
force was therefore exactly l/r^ beyond h) . Very roughly 
this corresponds to a P lumm er law smoothing e ~ /i/3 
(e.g. Springel & White 1999), although a Plummer law 
gives 1% force accuracy only beyond lOe. 

We have performed numerous tests of the code, among 
them tests of self-similar evolution of power-law spectra 
in critical density models and stable evolution of known 
halo profiles. The simulations took ~ 200 time steps from 
z = 60 to z = 3. Comparison of final particle positions 
suggested the time step criterion was conservative. We 
have additionally compared the TreePM code with a cos- 
mological Tree code (Springel & White 1999) and found 
good agreement in the clustering statistics for several dif- 
ferent initial conditions, including one of those used here 
(V. Springel,, private communication). 

With both the PM and TreePM codes we ran 3 real- 
izations of 4 models: the 'fiducial' ACDM power spec- 
trum, and 3 filtered versions with ko — 10/iMpc^^, closel y 
approximating Fig. 1 of Kamionkowski & Liddle (1999), 



ko = 5hMpc~^ and fco — 2/iMpc~^ which show a larger 
effect more easily resolved by these relatively small sim- 
ulations. For each of the 3 realizations the same ran- 
dom phases were used for all 4 power spectra to allow 
inter-comparison. As additional checks on finite volume 
and resolution effects we also ran simulations in boxes 
of side 50h~^ Mpc and 35ft.^^Mpc finding excellent agree- 
ment where the simulations overlapped. 



4. RESULTS 



4.1. 



Visual impression 

In Fig. |l] we show slices through the particle distribu- 
tions of our 4 models. The most extreme model, with 
ko = 2 h~^Mpc looks markedly different from the oth- 
ers, with smooth low density regions (the simulation ini- 
tial grid is still clearly visible), and a lack of substructure 
in the higher density areas. The differences between the 
other panels are more subtle, and in all cases, are really 
only apparent on the smallest scales. 

4.2. Power spectrum 

Most probes of small-scale structure, other than the ob- 
ject abundance, depend upon the non-linear power spec- 
trum. The process of gravitational collapse transfers power 
from large-scales to small, and can generate a k~^ tail 
in P{k) if it is absent initially. Fitting formulae for the 
non-lin ear p ower spectrum such as that of Peacock & 
Dodds ( |1996| ) are not applicable for spectra, such as ours, 
which have regions with n < —3. We use our PM and 
TreePM simulations to study the non-linear power spec- 
trum. 

As can be seen in Fig. ^ the scales of interest are non- 
linear by z = 3 and small-scale power removed by filtering 
has been regenerated by collapse of large-scale modes. The 



Slices though the particle distribution in the TrccPAI simulations of our 4 models at z = 3. From left to right, 
s show the models with filter scale ko = 2 h Mpc~^, ko = b h Mpc~^ (top row), fco = 10 /i Mpc~^, and fiducial 
(bottom row). The box side-length in each case is 25 /i~^Mpc and the slice thickness 0.25 /i~^Mpc. 
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Fig. 2. — The linear and non-linear power spectrum at 
z = 3. We show the mean from 3 runs of our TreePM sim- 
ulations for each of 4 models: our fiducial ACDM model 
(triangles), one filtered at ko = 10/iMpc~^ (open squares), 
fco = 5/iMpc~ (open circles), and ko = 2/iMpc^ (three 
pointed stars). The so lid line shows the prediction of Pea- 
cock & Dodds (1996) for the fiducial model (the theory 
is not applicable to the filtered models) and the dotted 
lines are the linear theory input spectra. While the N- 
body simulations have the same random phases and so are 
directly comparable, comparison with the analytic models 
requires error bars. We show the la error on the mean of 
the fiducial model computed from our 3 realizations. The 
errors on the other models are similar. 



fiducial model has good ag reeme nt with the fitting func- 
tion of Peacock & Dodds ( 1996| ) on intermediate scales, 



though for scales smaller than k ~ 10/iMpc~ in the 
TreePM simulations we obtain more power than Peacock 
& Dodds predict by a factor of about 2, independent of 
the realization or box size. We believe this is due to the 
very fiat nature of the linear theo ry spectrum on these 
scales: Jain & Bertschinger ( 1998 ) found a similar dis- 
crepancy with Peacock & Dodds for n = —2 spectra (see 
their Fig. 7). We note that other fitting formulae have 
been developed which it is possible would better repro- 
duce the be haviou r of our fiducial model (see e.g. Jain, 
Mo & White |1995| ; Ma |1998D . 

To focus on the dynamics of the power regeneration 
we show the evolution of the mass power spectrum for 
our fiducial model and one filtered model (with fcp = 
5/iMpc-^) in Fig. |. We use the average of 3 realiza- 
tions of PM simulation output here since the greater par- 
ticle density allows us to probe smaller amplitude fluc- 
tuations before shot-noise contamination becomes severe. 
The PM and TreePM simulations agree on the power up 
to fc ^ 20h Mpc~^ suggesting we resolve the relevant scales 
with our PM code. 

Notice that even at z = 6, the "peak" in power in- 
troduced by Eq. ^ has disappeared and small-scale power 
has been regenerated. The difference between the fiducial 
and filtered models grows progressively smaller as the evo- 
lution proceeds. For comparison, the generation of non- 
linear power has also been studied in numerical experi- 
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Fig. 3. — The non-linear power spectra as a function of 
redshift for our fiducial model (solid symbols) and one fil- 
tered model with ko = 5/iMpc~^ (open symbols). The 
symbols represent the mean of 3 runs of our PM simula- 
tions, the error bars show the la error on the mean esti- 
mated from our 3 realizations at z = 0. (However recall 
that for each model the simulations have the same ran- 
dom phases.) The spectra are scaled by (4a) ~^ to reduce 
the effect of linear evolution and highlight the non-linear 
growth. 
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ments by Little, Weinberg & Park ( 1991 ), who s t udied 
scale- invariant models, Melott & Shandarin (1990, 1991, 
1993 ) and by Bagla & Padmanabhan (1997), amongst oth- 
ers. 

Finally it is of interest to ask how the redshift space 
power spectra evolve. Typically the redshift space spectra 
appear closer to the linear theory power spectrum than 
the real space spectra. In Fig. |j we show the redshift 
space mass power spectrum as a function of redshift, as in 
Fig. H for the real-space spectra. We can see that even the 
redshift space spectra have a tail of power at small scales, 
induced by the non-linear clustering. 

4.3. Ly-a forest 

There has been a great deal of progress in theoretical un- 
derstanding of the Ly-a forest recently, due in large part 
to hydrodynamic simul ations (Cen et al. 1994; Zhang, An- 
ninos & Norman 1995: Miralda-Escude et al. 1996; Hern- 



Theuns et al 
1999|) . In 

at high z (> 



quist et al. 1996 ; Wadsley 

[19984 ^998b 
al ■ " ' 



Bondh99€; Z hang e t al. 1997 



Dave et al. [1998| ; Bryan et 
these simulations, it has been found that 
2), most of the absorption in Ly-a forest 
spectra is due to a continuous, fluctuating photoionized 
medium. The physical processes governing t his ab sorb- 
ing gas ar e simp le (see e.g., Bi and Davidsen 1997, Hui 
& Gnedin 



1997) 



and as a result, the optical depth for 
absorption at a particular point can be rel ated directly to 
the underlying matter density (Croft et al. 1997). Because 
of this, observations of the Ly-a forest in quasar spectra 
can be potentially very usefu l for probin g the clustering 



of matter (e.g Gnedin 1998, Croft et al. 1998, Nusser & 



Haehnelt gOOOj). 

We have generated simulated Ly-a forest spectra from 
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k (h/Mpc) 



our PM A^-body outputs at z = 3 in order to test how 
constraining Ly-a measurements could be for the mod- 
els described in this paper. To do this, we follow asim- 
ilar procedure to that outlined in Hui & Gnedin (1998) 
and Croft et al. (1998). We bin the particle distribution 
onto 512"^ density and velocity grids using a cloud- in-cell 
scheme, and smooth with a Gaussian filter of width one 
grid cell. We convert the density in each cell to an op- 
tical depth for neutral hydrogen absorption by assuming 
the 'Fluctuating Gunn-Pete rson Approximation" (FGPA) 
(see Croft et al. [ |1997| , |1998| ; see also Hui & Gnedin |1998[ ) 
and assign a temperature to the cell using a power-law 
density-temperature relation. In all our tests, we us e the 
form T = Top'-'-\ with 7 = 1.5 (see Hui & Gnedin |1997 



for the expected dependence of 7 on reionization epoch). 
We set the coefficient of proportionality between density 
and optical depth by requiring that the mean transmitted 
flux, (F) — 0.684, in a ccordance with the observations of 
McDonald et al. (199£). We run 256 lines of sight paral- 
lel to each of the three axes through one simulation box 
and create mock spectra from a convolution of the optical 
depths, peculiar velocities, and thermal broadening. The 
conversion to kms~^ from /i^^Mpc at z = 3 in this model 
is a factor of 112. 
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Fig. 4. — The redshift space non-linear power spectra as 
a function of redshift for our fiducial model (solid sym- 
bols) and one filtered model with fco = 5/iMpc~^ (open 
symbols). The symbols represent the mean of 3 runs of 
our PM simulations. The spectra are scaled by (4a) 
to reduce the effect of linear evolution and highlight the 
non-linear growth. 



Fig. 5. — (top): The ID power spectrum of the flux mea- 
sured from our simulations of the 4 models (lines), all with 
Tq = lO^'K. The observational results, also at 2; = 3, of Mc- 
Donald et al. (1999) are shown as points, (bottom): The 
ID power spectrum of the flux for the fiducial model with 
3 different values of Tq. 
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On small scales, the finite pressure of the gas will in 
detail modif y its c lustering (see e.g., Hui & Gnedin |l998 
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Fig. 6. — The dimensionless 3D power spectrum of the 
flux measured from our Ly-a PM simulations of 2 of our 
models: fiducial ACDM (filled triangles), and the filtered 
model with ko = 2/iMpc~^ (open triangles). The linear 
theory mass power spectra, scaled down by an arbitrary 
factor to match the flux power spectra, are shown as dotted 
lines. 
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Fig. 7. — Panel (a): The probability distribution function 
of the flux, P{F), measured from our Ly-a PM simulations 
of the 4 models, aU with Tq = lO^K. Panel (b): P{F) for 
the fiducial ACDM model, with three different values for 
the temperature at the mean density. 



Bryan et al. 1999| , Theuns et al. 1999| ), tending to make 



the gas density field smoother than the dark matter only 
outputs of our simulations. We have also imple mente d 
a 2-species version of "Hydro-PM" (Hui & Gnedin |1998| ), 
which takes these effects into account, but find that for our 
purposes here the main results are adequately reproduced 
by the pure PM runs. We expect that the details of the 
spectra that we create will also depend on other assump- 
tions about e.g., the reionization epoch and our simulation 
methodology. To the extent that we are interested primar- 
ily in relative comparisons between models this should not 
be cause for concern. 

For each set of mock spectra we compute the one- 
dimensional flux power spectrum, using an FFT, and show 
the results in Fig. ||. In the top panel of this figure, we have 
set the temperature of the gas at the mean density, Tq, to 
be equal to 10*K for all models. The different curves show 
the effects of linear power suppression at different values 
of ko, while th e poin ts are the observational results of Mc- 
Donald et al. ( |1999D . While a suppression in fiux power is 
seen in Fig. it is very small. If we vary the value of Tq, 
the change in the thermal broadening scale causes a more 
dramatic effect. This can be seen in the lower panel of 
Fig. ^ where we show results for the fiducial linear power 
spectrum only. In general, the flux power spectrum shape 
on small scales will depend on the temperature of the gas 
through thermal broadening and flnite gas pressure, as well 
as the non-linearity of matter clustering. In the context 
of this ACDM model, it seems as though the observations 
of McDonald et al. are consistent with a fairly high mean 
gas temperature, although a more detailed study involv- 
ing hydrodynamic simulations is needed to give definitive 
results. What is certain from the present study is that 
the one-dimensional flux power spectrum provides little 
constraint on our models with suppressed linear power. 
Clustering in the flux has apparently been regenerated by 
non-linear gravitational evolution in a similar fashion to 
that seen in Fig. There may be a good side to this, 
though, as insensitivity to the amount of small-scale lin- 
ear power will mean that estimates of the temperature of 
the IGM made by looking at small-scale clustering of the 
flux (as in Fig. |^) should be more robust than expected. 

If we assume isotropy of clustering, the three-dimensional 
flux power spectrum, A'j^{k) can be simp ly rec overed from 
the one-dimensional one (see C roft et al. 1998 for details). 
It was found by Croft et al. (1998) that on sufficiently 



large scales the shape of A|, measured from simulated 
spectra matches well that of the linear theory mass power 
spectrum, A^(fc). In Fig. || we test this using the model 
with ko — 2 h~^Mpc, and the fiducial model (both with 
To = lO^'K). We can see that there is not much difference 
between A|,(fc) for the two (the same is true of the two in- 
termediate models, which we do not plot). The linear the- 
ory mass power spectrum (arbitrarily normalized) is shown 
for comparison. On scales approaching the box size, cos- 
mic variance is large enough to account for the difference 
between the linear curve and the points. On smaller scales, 
there is still scatter, but the points taken from the simula- 
tion with less linear power are systematically a bit lower. 
We might expect the simulation points to start to trace 
the linear theory shape around the scale of non-linearity, 
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where A^(fc) becomes comparable to 1, which from Fig. ^ 
is around k ^ 1 — 2 h Mpc^^. If we look at Fig. this 
does seem reasonable, and we do find similar results even 
if we assume different gas temperatures. On smaller scales 
though, A|,(fc) has been regenerated by non-linearity, so 
that the exact relationship between A|,(fc) and the mat- 
ter clustering is complex, and as in Fig. ^ the differences 
between models small. 

Another statistic which we can check, in order to see 
whether suppression of linear power has caused any changes 
in higher-order clustering, is the probability distribution 
of the flux. We plot this in Fig. 0, showing the 4 mod- 
els with different linear power (and all with To = 10"' K) 
in the top panel. There are small differences between the 
models, particularly at the high flux end, where the mod- 
els with more power appear to have more truly empty 
regions. These small differences are likely to remain unob- 
served though due to the difficulty of accurate continuum 
fitting. This has implications for studies which use the flux 
PDF information to constrain the amount of line ar power 
on the Jean's scale (e.g., Nusser & Haehnelt 2000 ). For the 
same reason that the flux power spectrum does not change 
much on small scales (generation of power) , these methods 
are likely to also be insensitive to power truncation of the 
type we are considering. 

The lower panel of this figure shows results for our fidu- 
cial model with different gas temperatures. There appears 
to be little difference between the curves, although we have 
found that differences do appear if the spectra are sub- 
jected to a moderate amount of smoothing (e.g. with a 
50 km s~^ Gaussian; not shown). 

From our tests with both sets of statistics, we find that 
the Ly-a forest is not a promising discriminator between 
the models we are considering here. Two effects conspire 
to mask any differences in the Ly-a measurements on the 
small scales where there are large differences in the linear 
power spectra. First, the thermal broadening has the effect 
of smoothing the spectra (the thermal width of features is 
a few 10s of km s""'^). Second, non-linear evolution of the 
density field causes power to be rapidly transferred from 
large to small scales. For these models the scale of non- 
linearity at z = 3 at is about the same or larger than 
the scale at which there are large differences in the linear 
power spectra. On smaller scales, the shape of the three- 
dimensional flux power spectrum no longer follows that of 
the linear mass power spectrum. 

4.4. Halo abundance 

Ideally we would like to evolve a large volume and study 
the number density of small halos present today within 
a larger halo such as the Milky Way. This is not pos- 
sible with the limited dynamic range of the simulations 
presented here. While many effects could potentially dis- 
turb all the small halos as they interact with each other 
and a larger halo, very high resolution numerical simula- 
tions suggest that this may not be the case in practice. 
Moore et al. ( |l999a| ) find that a large number of small ha- 
los are not disrupted, so that the number remaining will 
still be a substantial fraction of the number that existed 
in the proto-galaxy. In this paper, we focus on the num- 
ber density of halos in our simulations at z = 3 where our 
25h~^M-pc box is just about to go non-linear. We assume 
that these small halos would become incorporated into a 
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Fig. 8. — The mass functions of our simulations, us- 



ing the HOP algorithm of Eiscnstein & Hut ( |1998D to 
find halos. We show results for fiducial model (trian- 
gles) and our 3 models with reduced small-scale power: 
fco = 10/iMpc^^ (open squares), fco = 5/iMpc~^ (open cir- 
cles), and fco = 2/iMpc~ (three pointed stars). We have 
co-added the mass functions of the three realizations to 
construct an "average" mass function - the error on the 
mean as calculated from the 3 realizations is smaller than 
the size of the symbols. The solid lines show the predic- 
tions of the Press-Schechter theory. We have used top-hat 
smoothing and Sc = 1-5 in the calculation of the Press- 
Schechter predictions. 
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Fig. 9. — The simulation mass functions, as for Fig. ^ 
except that a friends-of-friends scheme (with linking length 
e = 0.2) was used to find halos. 
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larger halo at later times by the usual evolution of cluster- 
ing, and that the fraction that survive disruption can be 
predicted by the re ferring to the detailed calculations of 
Moore et al. ( 1999a ). Here we are only interested in deficit 
in the number of small halos in our suppressed models rel- 
ative to the fiducial model. This relative fraction should 
be similar at the redshift of our simulation box to what it 
would be at 2 = 0, although the absolute number of halos 
could only be quanitified using simulations like those of 
Moore et al. Similar assum ptions to ours were also used 
by Kamionkowski & Liddle (199£). 

We show in Fig. H the Press-Schechter predictions for 
our fiducial and filtered models and the numerically de- 
termined mass functions. As there is no perfect algorith- 
mic definition of a "group" of points, the mass function is 
sensitive to a small degree to the halo finding algorithm. 
We h ave u sed both the Friends-of- Friends (FOF; Davis 
et al. 1985) algorithm, with l inking length 0.2, and the 
HOP (Eisenstein & Hut 1998) halo finding algorithm to 
construct these mass functions. We find that the mass 
function differed slightly if we changed the parameters in 
the algorithms or the algorithm used, and show results for 
both of these schemes. However, these differences should 
not affect our main conclusions, as we are interested in 
comparing models with different amounts of small-scale 
power to each other. 

The N-body mass functions and the Press-Schechter pre- 
dictions are shown in Figs. H and ^ for the HOP and FOF 
algorithms respectively. The "shelf" at the low mass end of 
the N-body mass functions arises because of the minimum 
number of particles allowed to form a group. There are no 
very low mass halos in the simulation. If we use Sc = 1-69 
with a top-hat window in the Press-Schechter predictions 
we find that the mass functions have too few large-mass 
halos compared to the HOP N-body results for both z = 4 
and 3. The two can be brought into better agreement if 
we decrease 5c to 1.5 as we have done. With this modi- 
fication, the Press-Schechter predictions overestimate the 
FOF results by a factor of up to 2. To check for simula- 
tion artifacts we also ran several larger boxes. We find that 
the mass functions from a sequence of larger boxes (up to 
50/i~^Mpc) with different random phases match smoothly 
and stably onto the mass function of these simulations, 
suggesting that there are no finite volume or sample vari- 
ance effects operating. As a final check, a completely sepa- 
rate analysis chain using a different N-body code (a cosmo- 
logical Tree code) obtains the same mass function at z = 3 
for one of our runs (V. Springel; private communication). 

Figs. H and ^ suggest that the number of halos is in- 
deed governed by the linear theory power spectrum. The 
amount of suppression relative to the fiducial model is ro- 
bust to the parameters of our group finding algorithm or 
the algorithm used. The absolute number of halos can 
in principle be predicted from statistics of the initial den- 
sity field, although there are uncertainties related to the 
definition of halos in the simulations and parameters in 
Press-Schechter theory. 

We find that in order to reduce the number of small 
halo s by a large factor (for example Kamionkowski & Lid- 
dle ( 1999D recommend about an order of magnitude), we 
require a fairly severe filtering of the fiducial model, using 
a filter with ko — 2 h Mpc~^. 

Finally we remark that this set of simulations does not 



have enough mass resolution to probe the structur e of the 
halos we find. However simulations by Moore et al. (1999b) 
suggest that the halo structure will not be sensitive to the 
filtering of the initial power spectrum. This lends some 
support to our assumption that the amount of disruption 
of the small halos when they become incorporated into a 
larger halo does not depend on the alterations we have 
made to the initial power spectrum. 

5. CONCLUSIONS 

While the essential picture of hierarchical formation of 
large-scale structure in a universe containing primarily 
cold dark matter appears to work well, some puzzles re- 
main. One of these is the paucity of dwarf galaxies in the 
local neighborhood. One resolution of this "lack-of-small- 
halos problem" is a modification of the initial power spec- 
trum, reducing the amount of small-scale power. There ex- 
ist inflationary models which can accomplish this, though 
the scale of the modification must be put in by hand. 
Other approaches, such as assuming that the universe is 
dominated by Warm Dark Matter (WDM) , will have a sim- 
ilar effect (and both approaches ma y solv e other problems: 
see e.g. Sommer-Larsen & Dolgov 19991 ). In models with 
reduced small-scale power structure forms in a top-down 
manner over a range of scales near the break, so ansatze 
developed for the "traditional" bottom-up scenario should 
be treated with caution. In this work we have used numer- 
ical simulations to address the question of how one could 
constrain such a modification of the initial power spec- 
trum. We note that we have dealt in detail only with a 
model with suppressed initial power. In a WDM model 
the deficit of power arises from the dark matter velocity 
dispersion, and so such as model may behave slightly dif- 
ferently, at least on the smallest scales. 

We find that the halo mass function depends primar- 
ily on the linear theory power spectrum, so a suppression 
of small-scale power does reduce the number of low mass 
halos. While the Press-Schechter theory predicts qualita- 
tively the right behavior, its free parameter (Sc) must be 
adjusted to fit the N-body results. To reduce the num- 
ber of 10^'~'Mq halos by a factor of > 5 compared to our 
fiducial model requires a fairly extreme filtering of the pri- 
mordial power spectrum, and the structure that forms in 
such a model appears qualitatively different to the fiducial 
ACDM model (Fig. [|). 

Collapse of large-scale structures as they go non-linear 
regenerates a "tail" in P{k) if it is suppressed in the initial 
conditions (and this holds in redshift as well as real space). 
Thus probes which measure primarily the evolved power 
spectrum are less sensitive to reduced small-scale power 
than one might think. We particularly examine measure- 
ments of clustering from the Ly-a forest flux. On the scales 
which govern the number of small halos, choosing a differ- 
ent gas temperature affects Ly-a clustering much more 
strongly than suppressing the linear power spectrum. The 
matter power spectrum measurement made from the low 
resolution Ly-a forest spectra by Croft et al. (1998) probes 
scales just above this, which are still linear, and offers es- 
sentially little constraint on these models. Any extension 
of these simple Ly-a forest measurements to smaller scales 
must necessarily have less general conclusions drawn from 
them. 

Given that the number density of collapsed objects 
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seems to be the most sensitive probe of this smaU-scale 
modification of the power spectrum, other observations 
which depend on this should be used to make consistency 
checks. At the moment, the obvious choices, such as the 
number density of damped Ly-a systems, or the redshift 
of reionization induced by the formation of the first stars 
and quasars, are difRcuh to predict accurately from the- 
ory. Their potentially strong discriminatory power will 
make them useful eventually though, as we learn whether 
more of Cosmology's puzzles can be resolved by an absence 
of sm all-scale power. 
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